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Abstract. We theoretically explore the optical flux lattices produced for ultra- 
cold atoms subject to laser fields where both the atom-light coupling and the 
effective detuning are spatially periodic. We analyze the geometric vector 
potential and the magnetic flux it generates, as well as the accompanying 
geometric scalar potential. We show how to understand the gauge-dependent 
Aharonov-Bohm singularities in the vector potential, and calculate the continuous 
magnetic flux through the elementary cell in terms of these singularities. The 
analysis is illustrated with a square optical flux lattice. We conclude with an 
explicit laser configuration yielding such a lattice using a set of five properly 
chosen beams with two countcrpropagating pairs (one along the x axes and the 
other y axes), together with a single beam along the z axis. We show that this 
lattice is not phase-stable, and identify the one phase-difference that affects the 
magnetic flux. Thus armed with realistic laser setup, we directly compute the 
Chern number of the lowest Bloch band to identify the region where the non- 
zero magnetic flux produces a topologically non-trivial band structure. 



PACS numbers: 67.85.-d, 32.10.Fn, 33.60.+q, 37.10.Gh, 67.85.Hj, 37.10.Vz 
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1. INTRODUCTION 

Atomic quantum gases are systems where condensed matter and the atomic physics 
meet. Cold atomic gases exhibit a number condensed matter phenomena [TJ [21 
|3J |3], such as the superfluid-Mott transition [5], Berezinskii-Kosterlitz-Thouless 
superfluidity 6], and the Bose-Einstein condensation to Bardeen-Cooper-Schricffcr 
(BCS) crossover [3 |8]. Because the atoms comprising these quantum gases are 
electrically neutral, no vector potentials affect their center of mass motion. Such a 
vector potentials might provide the Lorentz force essential for the magnetic phenomena 
in solids, such as the quantum Hall effect }9j. The standard way to produce an artificial 
magnetic field is to rotate an atomic cloud leading to a non-trivial vector potential 
in the rotating frame of reference (TU1 E] ■ The various proposed schemes to create 
effective magnetic field for ultra-cold atoms without rotation [12] can be divided into 
two categories. 

The first relies on a primary optical lattice which traps atoms at its sites. 
The magnetic flux is created by inducing asymmetric tunneling between lattice 
sites, so that atoms acquire a non-zero phase after completing a closed loop along 
a plaquette [HI HH [HI HSJ [17l [HI Ell [20l EU [22]. Such asymmetries can be 
induced by laser-assisted tunneling [HI [HI [HI (HI [HI EU [23] or using time-dependent 
lattices [HI H3 HH HOI 122] • 

The second group of proposals is based on the concept of geometric gauge 
potentials which are encountered in many areas of physics [231 1221 [23 123 HH [23 
l3"01 I3T1 l3"2"j . In atomic gases, the geometric vector and scalar potentials were first 
considered in the late 90's for atoms interacting with the laser fields [331 1331 l35j . 
where the atoms are "dressed" by laser beams. The resulting position dependence of 
the dressed internal states leads to geometric vector and scalar potentials. The method 
can provide a non-zero effective magnetic field using non-trivial spatial arrangements 
of laser fields [MlETlEHlESSniEIlSSES] or position-dependent detuning of the 
atom- light coupling [331 [45j [46] . In these approaches the magnetic flux through the 
atomic cloud scales linearly with the cloud's extent [40l [42l [44j [45] , not it's area. For 
large systems this is a major obstacle in reaching the sizable magnetic fluxes required 
to achieve the fractional Hall effect [4"7] . 

A new class of geometric potentials termed "flux lattices" were recently shown 
to yield a magnetic flux proportional to the surface area of the atomic cloud [48] , see 
also |49j . In this proposal, a two-level atom was coupled to a spatially periodic laser 
field where both the atom-light coupling and the detuning term were oscillatory. This 
approach simultaneously generates a non-staggered magnetic flux along with a lattice 
potential thus providing an optical flux lattice [4"8] . 

The vector potential A plays an important role in the quantum physics |50] . It is 
featured in the Pearls substitution [5U [52] [53] widely used in the tight binding models 
in solids to describe the motion of charged particles in the magnetic field. Specifically, 
the tunneling matrix element between the lattice sites and rg acquires the Pearls 
phase factor proportional to J r ^ A • dr. Similar Pearls phase factors emerge also in 
the tunneling matrix elements between the sites of an optical lattice for electrically 
neutral ultracold atoms affected by the artificial magnetic field. 

It is instructive that the geometric vector potential for the optical flux lattices 
contains gauge-dependent singularities. To avoid these singularities, Ref. [3E] 
concentrated on the magnetic flux rather on the vector potential. Here we explore 
optical flux lattices explicitly in terms of their geometric vector potential (which 
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generates the magnetic flux) and an accompanying geometric scalar potential. We 
show how to understand the gauge-dependent Aharonov-Bohm (AB) singularities [5U] 
appearing in this vector potential and we calculate the continuous magnetic flux 
through the elementary cell in terms the singularities. 

Subsequently, we analyze a square optical flux lattice and describe a way of 
creating it using a combination of state-dependent lattice potentials and Raman 
transitions. A related setup proposed recently by Cooper and Dalibard aimed at 
producing triangular and hexagonal optical flux lattices used three coplanar lasers 
intersecting at 120 degrees, with an additional beam normal to the plane spanned by 
the first three [54] . Here we present an explicit laser configuration yielding a square 
flux lattice and directly compute the Chern number of the lowest Bloch band. We 
identify the region where the non-zero magnetic flux produces a topologically non- 
trivial band structure for this lattice. 

2. Hamiltonian and its eigenstates 

Before focusing on a specific physical system, we begin by considering the very general 
problem of a multi-level atom moving in the presence of a spatially inhomogenous 
coupling Hamiltonian (for example produced by a combination of optical and magnetic 
fields). The Hamiltonian describing such a combined internal and center of mass 
motion is 

i+M(r), (1) 

where p = — i/iV is the atomic momentum; 1 is the identity operator; U (r) is a state- 
independent "scalar" potential; and M is the state-dependent part of the Hamiltonian. 
Here we focus on the case where the atom affected by the light fields behaves like a 
spin in a magnetic field, so the state dependent Hamiltonian M(r) is 

M(r) = fl F = OFh, (2) 

where the vector Jl = fi(r) = (f2 a; (r) ) Cl y (r), O z (r)) describes the spatially dependent 
coupling between the atomic internal states, f2(r) = |S7(r) j being the total coupling 

strength; F = (^F x ,F y ,F z ^j is a vector operator satisfying the angular momentum 

algebra (to be referred to as the spin operator); and Fq = -Fn(r) is the spatially 
dependent projection of F along fi. The position dependence of the Hamiltonian 
AT(r) therefore originates from the position-dependent atom-light interaction through 
the coupling vector f2(r): a rapidly varying effective magnetic field. The physical 
implementation of such a Hamiltonian - physically equivalent to the Zeeman effect for 
a spatially dependent magnetic field - will be discussed in Sec. [5] 

Figure [lj, depicts the (quasi-)spin-l/2 case, where HQ Z is the light-induced 
detuning between the two internal atomic states \mp = ±1/2), and H (£l x ±iQ y ) /2 
is the transition matrix element coupling the two states together. In what 
follows, we do not restrict ourselves to the spin-1/2 case: for N internal states 
{\t71f} \mp = —f, f + 1, ... /}, the quantity / = (N — l)/2 is the total angular 
momentum quantum number. The spin-1 case depicted in Fig. [TJd, has / = 1 and 
N = 3. We emphasize that the atomic states |m.F) do not necessarily represent the 
true spin states. They can be the atomic internal states of arbitrary origin, provided 
that the operator F featured in the atomic Hamiltonian obeys the angular momentum 
algebra. 



H 



P 

2m 



■U(t) 
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Figure 1. Traditional representation of M = f2 • F coupling scheme for (a) 
total angular momentum / = 1/2 and (b) total angular momentum / = 1. The 
bare atomic states are labeled by \mp), and in both cases are detuned from each 
other by a frequency Q z . These levels are coupled with strength proportional to 



2.1. Diagonalization via a unitary transformation 

The projected momentum operator Fq entering the coupling Hamiltonian M can be 
related to the eigenstates of F z via a unitary transformation Ssi where 



with 



Fn — SqFzSq 1 



§ _ e -iF z <l>/h e -iF v 9/n e iF z 4,/h_ 



(3) 



(4) 



We parametrize the coupling vector fi = (£l Xl fl y , f2 z ) in terms of the spherical angles 



tan 6 — 



f2y 



and 



cos# = 



0* 

n 



(5) 



shown in Fig. [2] 

The operators F 2 and F z have eigenstates \m F ) = \ f, mp) identified by the total 
angular momentum / and its e z projection m F with eigenvalues 

F 2 \m F ) — h 2 f(f + 1) \rn F ) 
F z \m F ) = Hm F \m F ) . 



Multiplying the last equation by Sq, one has 

Fq \ fh F ) — hm F \ fh F ) , where 



Imp) = Sn Imp) 



(0) 



Therefore the coupling Hamiltonian M = £IFq has a set of position-dependent 
eigenstates \rhp) = \rhp(r)} related to the eigenstates of F z via the position-dependent 
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Figure 2. Representation of the coupling vector fl = (Q x , (l y , fl z ) in terms of 
the spherical angles 9 and <f> . 



unitary transformation Sq = 5n(r). Using Eq. Q for So and the fact that \m F ) is 
an eigenstate of .F z , one arrives at the eigenstates 




and energies 

v mp = hm F n. (8) 

of the coupling Hamiltonian M. Here, V(r) = V mF is the position dependent energy 
of the local eigenstate |mF(r))- Interestingly, similar kinds of eigenstates give rise to 
artificial gauge potential terms describing the rotation of diatomic molecules [26l EI] 
and the physics of atomic collisions [SSj [S5] . 

3. Gauge potentials 

For an atom subject to the Hamiltonian of Eq. 0, the state- vector describing both 
it's internal and motional degrees of freedom can be expressed in the basis of dressed 
states 

l*(r,t)) =J2^m' F (r,t)\m' F ), 

m' F 

where ip m ' F ( r , t) is a wave-function describing the atom's motion in the basis of local 
eigenstates \m' F ) = \rh' F {r)). We are interested in the situation where f2 > 0, so the 
local eigenstates \rh' F ) are everywhere non-degenerate. 

If an atom is prepared in one of these dressed states with m' F = m F , and 
it's characteristic kinetic energy is small compared to the AE — Ml the energy 
difference between adjacent spin states, the internal state of the atom will adiabatically 
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follow the dressed state \m F ) as the atom moves, and contributions due to other 
states with m' F ^ m F can be neglected. Projecting the full Schrodinger equation 
ih t)) = H\ty(r,t)) onto the selected internal eigenstate \m F ) yields a reduced 
Schrodinger equation for the atomic center of mass motion ihip mF (r, t) — Hip„ lF (r, t) 
with an effective Hamiltonian 

H= [p -A(r)} + + + 
2m 

In this expression, the geometric vector 

A=A(r) =ih(rh F \ V\m F ) (10) 

and scalar 

W = W(t) = ^. J2 l(^|V|m F )| 2 (11) 
2m ^-^ 

potentials appear due to the position dependence of the atomic dressed states. The 
vector potential can be interpreted as the average center of mass momentum of the 
selected internal state \mp) = \mp (r)). The scalar potential W(r) emerges due to the 
elimination of the remaining atomic internal states. It represents the kinetic energy 
of the oscillatory micromotion [57J [55] due to the tiny transitions to the eliminated 
states \rh F ) — \rh F (r)) with m" F ^ rhp. 

Equation ^ contains three distinct scalar potentials: (a) the state independent 
potential U(r) featured in the initial Hamiltonian (JlJ, which we shall call U(r) the 
"scalar light shift;" (b) the "adiabatic scalar potential" V(r) = V mF arising from 
spatial variations in the magnitude O(r); and (c) the "geometric scalar potential" 
W(r) described above. All three contribute to the potential energy of atoms in the 
dressed state basis. 

3.1. Vector and scalar potentials 

Using Eqs. ([8]), the matrix elements featured in the vector and scalar potentials are 
(m' F \ V \rh F ) = l - e <^~^)4> ( m » | ]^ mp _ F ^ v</) _ p y V9~\ \m F ) , (12) 

with F z = exp(iF y 8/h)F z exp(—iF y 6/h) — F z cos 6 + F x sin 9. Using the identities 



(mp\F x \mF) — (mp\ F y \m F ) — and (mp\F z \mp) = mp, Eqs. ( 10 )-( 12 1 provide 
the vector potential 

A{r) = hm F (cos 9-1) Vcj>. (13) 
The vector potential (|13[) is maximum in magnitude when mp = ±/, and is zero for 



mp = 0. For mp = 1/2, Eq. (13 1 reduces to the result presented in Ref. |12) . 

To determine the scalar potential, we need the off diagonal matrix elements of 
(m F \ V \ fhp ) which are 

(m' F \V\m F ) = -^e 1 ^-™'^ {m' F \ (P x sin 6V<f> + F y V^) \m F ) , (14) 
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for rh'p 7^ in p. Combining Eqs. (11) and (14), and using the completeness relation, 
we arrive at the geometric scalar potential 



W(r) = ^g ttVU 



with 



sm 2 e (V4>y + (yey 



(15) 



9f,m F = /(/+ 1) 



vp, 



(16) 



where we made use of: (ra^l F x \ijif) = (rnp\ Fy \mp) = gf iin , F /2. In particular, for 
/ = 1/2, one has g/, mF = 1/2, and the scalar potential reduces to that presented in 
Ref. [T2]. Generally W(r) depends both on the total spin / and also on its projection 



mp. For instance, for / = 1, one has gf,m F 
for rriF = and is a half of that for mf = ±l 



2 — mp, showing that W^(r) is maximum 



3.2. Alternative gauge 

Because each of the local eigenstates can be assigned an arbitrary position dependent 
phase \m' F ) = exp [i(fm F (r)] \rrip}, effecting a state-dependent gauge transformation, 



the vector potential «4.(r) in Eq. (13) is not unique. For example, when 



\m F ) 



(17) 



the initial vector potential given by Eq. ( 13 ) becomes 

A!(r) =hm F (cos 8 + 1) V<t>. 

This seemingly esoteric change can have significant impact because the 
contribution may be singular when cos 6* = ±1 (near the z axes, see Fig. [2]) if the 
factors (cos0 — 1) or (cos 8 + 1) do not compensate the singularity by simultaneously 
going to zero. The vector potential .4.(r), may have singularities when 8 = it, but 
in the alternative gauge the vector potential A!(r) has potential singularities when 
8 = 0: at spatially different points than the initial gauge! 



3.3. Magnetic flux 

The singularities in the vector potential correspond to AB type flux tubes (piercing the 
&x-&y plane) each with an integer flux quantum. Since the AB type flux containing an 
integer number flux quanta can not be observed [50j . the two vector potentials A.(v) 
and A!(y) are equivalent and produce the same effective magnetic field 



8(r) = V x A(r) =hm F V (cos 8) x Vq 



(18) 



The gauge-dependent AB singularities (if any) present in the vector potential must be 
absent in Eq. (18) for B(r). 

It is convenient to represent the magnetic flux density in terms of the unit vector 

n = n/n 

. VN X x VN y 
-hrriF 



S(r) 



N, 



(19) 



Thus if Q z alternates the sign at the points where tt x — Q y — 0, this might compensate 
the alternation of the sign of WN X x ViV y at these points, giving a non-zero magnetic 
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flux, such as the one given by Eq. (30 1 below. This shows the necessity to have an 



oscillating detuning f2 2 in addition to the oscillating coupling Q x + itt y 



As we show in Sec. 4.2 the geometric scalar potential W(r) contributes most 
significantly to the overall scalar potential U(r) + V(r) + W(r) at the maxima of the 
effective magnetic field where Q x + ifl y = 0, and is zero at the points of the minimum 
magnetic flux where Cl z = 0. 

3.4. Periodic atom-light coupling 

Given this general background, we now consider the case where the coupling vector 
ft = (Q x , ri y , fl z ) is spatially periodic in the e x -e y plane 

^(r + r„, m ) = fi(r), r n<m = na l + ma 2 , (20) 

where ai and a.2 are the primitive vectors defining a 2D lattice in the e x -e y plane, 
with {n,m} E Z. In this case, both the atomic internal dressed states |m^(r)) 
and the corresponding geometric potential .4.(r) have the same periodicity [usually 
the geometric scalar potential W(r) has a twice smaller periodicity than the initial 
Hamiltonian and the geometric vector potential] . Due to the periodicity of the vector 
potential the total flux over the elementary cell is zero 

' ' A dv= \ ff Stot ■ dS = 0, (21) 



^ /cell ^ J J cell 

where 2? to t = S(r) +Bab(i") is the total magnetic flux density with contributions both 
from the continuous (background) magnetic flux density B(r) and possibly a set of 
gauge-dependent singular fluxes of the AB type represented by Bab (r) . 

Thus it is strictly speaking impossible to produce a non-zero effective magnetic 
flux a over the elementary cell using the periodic atom-light coupling. However, 
this does not preclude a non-staggered continuous magnetic flux density 2?(r) over the 
elementary cell as long as the vector potential contains (gauge-dependent) singularities 
of the AB type carrying together a non-zero number of the Dirac flux quanta. The 
AB singularities are associated with the points where the Q x + iVL v goes to zero and 
hence cosO — > ±1. Deducting these non-measurable gauge dependent singularities, 
the remaining flux over the elementary cell can be non-zero 



Iff BdS = -l ff B AB (r)-dS. 

n JJcc\\ n JJccl\ 



(22) 



h 

The physical flux can hence be expressed in terms of the vector potential 

a' = -i^/ A-dv = -j^<l A'-dr, (23) 



singul /fc ^singul 

where the summation is over the singular points of the vector potential (emerging at 
cos —¥ — 1 for A and at cos 9 — !• 1 for A') around which the contour integration is 
carried out. In the neighborhood of each singular point (different for A and A) the 
vector potentials have a Dirac-string (the AB singularity) piercing the e x -e y plane, 
giving 

A^f -2hm F V<j), and A ->• 2hm F Vcj). (24) 
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Thus each integral in Eq. ( 23 1 provides an integer number of the Dirac flux quanta. 
To obtain a non-zero flux a' , the sum of the singular contributions must be non-zero. 
The flux is maximum if all these singular contributions have the same sign, as is the 
case for the square flux lattice considered below. 

To summarize, the optical flux lattice contains a background non-staggered 
magnetic field 3 plus an array of gauge-dependent Dirac-string fluxes of the opposite 
sign as compared to the background. The two types of fluxes compensate each other 
so the total magnetic flux over an elementary cell is zero as is required from the 
periodicity of the Hamiltonian. However, the Dirac-string fluxes are non-measurable 
and hence must be excluded from any physical consideration. As a result, a non- 
staggered magnetic flux over the optical flux lattice is possible. 

4. SQUARE OPTICAL FLUX LATTICE 

We now construct a simple model flux lattice generated by a spatially periodic coupling 
vector fi = (£l x ,Q y ,Q z ) with components 

Q :c = fl± cos(iE7r / a) , 

fly = fl± cos(j/7r/a), (25) 
Sl z = f2|| sin(x7r/a) sin(j/7r/a). 

This coupling has period 2a along e x and e y . For simplicity we define dimensionless 
coordinates x' — irx/a, y' — ny/a, and z' — irz/a. In the symmetric case, fl± = Oy, 
the scheme reduces to the one considered previously in Ref . [38] . As will be discussed 



in Sec. 5.1.2 the coupling vector in Eq. (25) can be produced for alkali atoms using 



five laser beams intersecting at right angles: a pair counterpropagating along e x , a 
second pair counterpropagating along e y , and a single beam propagating along e z . 



The total Rabi frequency resulting from Eq. (|25j is 

n 



+ (ni - nf) {fl + !l) + n\PJl (26) 

where /„ = cos(m'). The resulting adiabatic energies V mF = hmpQ have periodicity 
a, half that of the atom-light coupling. 

When S!^ > ftjj/2, the minima of the mj? < adiabatic scalar potential V(x, y) 
are positioned at x' n = nn and y' m = 7rm (brown dots in Fig. [3]) where Q z — 0. The 
energy maxima are positioned at x' nmax — ir(n + 1/2) and y' m max = 7r(m + /2) where 
the atom-light coupling vanishes: £l x + i£l y — > 0. Thus one has 

Emm — V(Trn,irm) — hrriF\/2Q,±, 

E m ax = V(n(n + 1/2), vr(m + /2)) = hmp^ . (27) 
In the vicinity of the energy maxima one has 

O^-fiL (X> - <, mQX ) (-1)", 

n y « -fi ± (y' - y' m . max ) (28) 

n z ps ft||(-i) n+m . 

Thus for odd (even) values of n + m the angle </> rotates clockwise (anti-clockwise) 
around the singularities of the vector potential positioned at x 1 = x' n max and 
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Figure 3. Sites of the square optical flux lattice corresponding to the minima 
of the adiabatic energy V(x, y) for mp < 0. 



v 

n- 



= y' nmax , whereas fi z alternates its sign when going from even to odd values of 
m. This ensures a non-zero magnetic flux over the elementary cell when integrating 



the vector potential around its singular points in Eqs. (23) and (24) 



4-1- Magnetic flux 

Consider now the flux passing through the elementary cell with x' G [0, 27r) 
y' G [0,2-71-). The vector potential A has Dirac-string singularities for Cl z = 
corresponding to odd values of n + m in Eq. ( 28 1 
two points are positioned at (n = 1, m 
magnetic flux quanta. In fact, integrating the vector potential around each singular 
point, Eqs. (23)-(24) yield the background magnetic flux 



and 

-fill 

Within the elementary cell these 
0) and (n — 0, m — 1), each providing Imp 



-y 



A - dr 



(29) 



singul 



In particular, for the spin-1/2 case {mp = 1/2) a measurable continuous flux over the 
elementary cell accommodates two Dirac quanta |48j . The same gauge independent 
magnetic flux a 1 is obtained using the alternate vector potential A 1 which contains 
gauge dependent AB singularities at different points: n = m = and n = m = 1, 
again each carrying Imp Dirac flux quanta. 



Using Eq. (19), one arrives at the explicit result for the magnetic flux density 



B(r) 



hmp ( — 
V a 



P {flf, 



3/2 



(30) 
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Figure 4. Rectangular coordinate flux lattice computed for /3 = 1. (a) Position 
dependent effective Zeeman magnetic field fl: the vectors denote the components 
in the e x -e y plane and the color depicts the e z component, (b) Geometric 
potentials. The color indicates the geometric scalar potential W(r), and the 
arrows denote the vector potential «4(r). (c) Magnitude of effective geometric 
magnetic field 13 along e z . Notice that regions of largest correspond to the 
maxima of W(r). 



where f u — cos(-u'), g u = sin(u'), and j3 — fly/fix < Equation ( 30 1 explicitly 
demonstrates that the magnetic flux, while non-uniform, is non-staggered, and its 
profile can be tailored by changing the ratio of the Rabi frequency amplitudes f3. 
It is evident that the magnetic flux is zero at the potential minima x n — na and 



y m — ma for finite values of /3. For (8 = 1, Eq. (30 1 is equivalent to a result obtained 
independently by J. Dalibard [55] , 

Figure [4] displays the spatial distribution of the effective Zeeman field fl, the 
geometric vector and scalar potentials A. and W, as well as the geometric magnetic 
field for (8 = 1. Figure [5] presents the geometric magnetic field 3 for various values of 
(9 showing that the most uniform magnetic field is reached for /3 = 1. 

4-2. Scalar potential 



Assuming flu = Q±, the second term entering the scalar potential [Eq. (IS)] is 

0(V0) -W hTfjm (B+wy (31) 



sin 2 ■ 



It is evident that sin 2 9 (V<^>) is zero at the minima of the adiabatic energy x n = na 
and y rn — ma. Additionally, it equals {it /a) at the maxima of the adiabatic energy 
x n ,max = na + a/2 and y m ,max = ma + a/2. This part of the scalar potential 
behaves similar to the adiabatic energy E_ (x,y), thus increasing the energy maxima 
by h 2 n 2 /8ma 2 . 



The first term entering the scalar potential [Eq. ( 15 )] is 



{V9) U (i + / 2 / 2 )(/ 2 +/ 2 ) ' (32) 

The gradient (V#) 2 is zero at the at the minima of the of the adiabatic energy where 
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Figure 5. Dependence of geometric gauge field B(r) on (5. (a) f) = 0.2 (b) 
P = 0.5 (c) /3 = 1.0 (d) p = 2.0 



9x = g y = 0, but is equal to (n/a) if f x = f y = 0, i.e. at the center of each plaquette 
thus raising the potential there. 

In this way, the geometric scalar potential is given by Eq. ( 15 ) together with 
Eqs. (31 1-(32 1. It is zero at the corners of a plaquette and reaches its maximum values 
at the center of the plaquette, thus behaving similar to the effective magnetic field, 
as is evident in Fig. [4] The scalar potential thus repels atoms from the area of high 
magnetic field at the center of the plaquette. 



5. Alkali atoms and light shifts 



In this Section, we shall first demonstrate a possible way to engineer the state- 
independent potential U(r) together with the state-dependent potential M featured 
in the general atomic Hamiltonian given by Eq. Subsequently we shall analyze 
atom-light configurations providing the square optical flux lattices. 

Let us consider a system of ultracold alkali atoms in their electronic ground state 
manifold illuminated by one or several laser fields which non-resonantly couple the 
ground states with the lowest electronic excited states. In the presence of an external 
magnetic field (yet without including the contributions due to the laser fields), the 
Hamiltonian for the atomic ground state manifold is 

H =H k + A M l3 + ^B-(gj3 + 9l l) . 
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where Hk — p 2 /2m is the kinetic contribution to the Hamiltonian; Ahj is the magnetic 
dipole hyperfine coefficient; and /is is the Bohr magneton. The Zeeman term includes 
separate contributions from J = L + S (the sum of the orbital L and electronic 
spin S angular momentum) and the nuclear angular momentum I, along with their 
respective Lande-g factors. We next consider the additional contributions to the 
atomic Hamiltonian resulting with off-resonant interaction with laser fields. 

As was observed in Refs. [60l EU [62] , conventional spin independent (scalar, U s ) 
optical potentials acquire additional spin-dependent terms near atomic resonance: the 
rank-1 (vector, U v ) and rank-2 tensor light shifts [50]. For the alkali atoms, adiabatic 
elimination of the excited states labeled by j = 1/2 (Dl) and j = 3/2 (D2) yields an 
effective atom-light coupling Hamiltonian for the ground state atoms (with j — 1/2): 



Hi 



The rank-2 term is negligible for the parameters of interest and hence is not included 
in Hl. Here E is the optical electric field; u v = — 2w s Aps/3(w — ujq) determines the 
vector light shift; Aps = ^3/2 — ^1/2 is the fine-structure splitting; hoJi/2 and 
are the Dl and D2 transition energies; and uio = (2wi/ 2 +W3/2)/3 is a suitable average. 
u s sets the scale of the light shift and proportional to the atoms ac polarizability. 

The contributions from the scalar and vector light shifts featured in Hl can 
be independently specified with informed choices of laser frequency oj and intensity. 
Evidently, the vector light shift is a contribution to the total Hamiltonian acting like 
an effective magnetic field 

ra„(E*xE) 

% = 

that acts on J and not the nuclear spin I. Instead of using the full Breit-Rabi 
equation [53] for the Zeeman energies, we assume that the Zeeman shifts are small in 
comparison with the hyperfine splitting - the linear, or anomalous, Zeeman regime - in 
which case, the effective Hamiltonian for a single manifold of total angular momentum 
F = J + 1 states is 

VB9F fT3 , „ ■, A , A h { 



H + H L = Us (E*-E) + ^ (B + B off )-F + -f (f 



Notice that B c g acts as a true magnetic field and adds vectorially with B, and since 
\gi/gj\ — 0.0005 in the alkali atoms, we safely neglected a contribution — fisgi^cft T/fi 
to the atomic Hamiltonian. We also introduced the hyperfine Lande g-factor 

/(/ + l)-^(^ + l)+J^0^ + l) 
9F gj 2/(/ + l) 

In 87 Rb's lowest energy manifold with / = 1, for which j = 1/2 and i = 3/2, we get 
gp = —gj/4 w —1/2. In the following, we always consider a single angular momentum 
manifold labeled by /, and select it's energy at zero field as the zero of energy. 

5.1. Bichromatic light field 

By combining state-dependent optical lattices along with "Raman coupling lattices," 
it is possible to create lattice potentials with large, non-staggered, artificial magnetic 
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(b) Level diagram 




Figure 6. Proposed experimental geometry, (a) Laser geometry for creating 
flux lattice showing the four circularly polarized beams in the e^-e^ plane 
with frequency u along with the linearly polarized beam traveling along e z 
with frequency ui + Suj. (b) Physical level diagram for three-level total angular 
momentum / = 1 case with mp states labeled, as is applicable for the common 
alkali atoms 7 Li, 23 Na, 39 K, 41 K, and 87 Rb. For reference, the diagram shows the 
decomposition of these optical fields into a± and 7r, but as discussed in the text, 
this is not an overly useful way of considering this problem. 



fields [48] even for alkali atoms [54j . Consider an ensemble of ultracold atoms subjected 
to a magnetic field B = B e z . The atoms are illuminated by several lasers with 
frequencies u and uj + Slo, where 5uj ~ |5_fMb^o/^| differs by a small detuning 
5 = gp^-BBs/h—Sta from the linear Zeeman shift between uif states (where \S\ <C Slu). 
In this case, the complex electric field E = E w _ exp(— iujt) + E w+ exp [— i(ui + 5ui)t] 
contributes to the combined magnetic field, giving 

B + B cff =B e z + [ (E* xE u _) + ( E* x E w+ ) 

+ (e:_xE u+ ) e~ i5ut + (e;xE w _) e iSut ] . 

The first two terms of B e ff add to the static bias field Boe z , and the remaining 
two time- dependent terms describe transitions between different wif levels. Provided 
Bq ^> |B e ff| and Sua are large compared to the kinetic energy scales, the Hamiltonian 
can be simplified by time-averaging to zero the time-dependent terms in the scalar 
light shift and making the rotating wave approximation (RWA) to eliminate the time- 
dependence of the coupling fields. The resulting contribution to the Hamiltonian 

H KWA = U(r)i + Sl-F, (33) 

takes the form of Eq. ([I]) once we identify the scalar potential 

U(r) = u s (e*_ -E w _ + E^ + -E w+ ) (34) 
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and the RWA effective magnetic field 



S + i 



h 



h 

fin 

-Re 



E* xE„ 



E* xE, 



e: + xe u+ 



(35) 



and ie y terms would 



This expression is valid for gp > (for gp < the sign of the e 
both be positive, owing to selecting the opposite complex terms in the RWA). The 
final form of this effective coupling shows that, while it is related to the initial vector 
light shifts, SI is composed of both static and resonant couplings in a way that goes 
beyond the restrictive B e g oc zE* x E form. This enables flux lattices in the alkali 
atoms. 

Importantly for practical flux-lattice configurations, ft z depends both on the static 
magnetic field and also the component of B e ff along e z . For practical considerations 
it is undesirable that the resonance condition be a function of the laser intensity, so 
we seek solutions without a contribution from this term. 



5.1.1. Two Raman beams First consider the straightforward example of the two 
counter propagating Raman beams used in existing experiments [H3II3IB4, 65, 6 61 IB"?]. 
In this simple case: 



E M = Ee ikRX e„, and E^ = Ee~ ikRX e z 



describing the electric field of two lasers counterpropagating along with equal 
intensities and crossed linear polarization, where — 2ir/\ is the single photon 
recoil wave vector, and Er = h 2 k R /2m is the associated recoil energy. The resulting 
scalar light shift U(r) and the effective magnetic field f2 describing the vector light 
shift are 



U(r) = u s (E* _ -E w _ + E^ + -E w+ J = 2u s E l 

11 = Se z + Or [sin (2knx) e x — cos (2kftx) e y ] , 

where FIr — u v E 2 /h. These describe a constant scalar light shift along with a 
spatially rotating effective magnetic field, as discussed in Ref. [64] which produced an 
artificial spin-orbit coupling, and in different notation, is equivalent to the proposal 
of Ref. [JD]. Because this Hamiltonian is only invariant under spatial translations 
with primitive vector u = Tr/kRe x , it would be expected to describe a periodic 
lattice. However, transforming the complete Hamiltonian according to the rotation 

U{x)HW (x), with U(x) = exp iF z (2kRX — tt/2) jh , completely removes the spatial 
periodicity. Instead, the transformed Hamiltonian becomes 

in which the position dependence has vanished from coupling vector Se z + Q,Re y . 
In this example, all spatial dependance (including the initial lattice structure) has 
been eliminated from the Hamiltonian in exchange for a matrix valued (though 
Abelian) gauge field. An additional spatially uniform radio-frequency magnetic field 
added to the mix forces the spatial structure to remain, creating a composite lattice 
potential (is . 
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5.1.2. Flux lattice configuration Next, we analyze the configuration shown in Fig. [6] 
where four beams with angular frequency uj intersecting at right angles in the e^-e^ 
plane, are joined by a fifth beam with angular frequency uj + Suj traveling along e z . 
The total electric field from these five beams are 

E w _ =E X + + E x - + E y + + E y - , and E w+ = E z 

where 

E x+ =E xy (V^ /2 cos0 p e z +e^ /2 sin^e y ) e «*-/a e i '*-»/ a e < * i,a 
E x - =£ XJ/ ( e -^ /2 cos(9 p e z - e^ /2 sin 6» p e y ) e -«0*/ V^^e" 4 ^ 
E y + =E xy (e^/ 2 cos9 p e z ~ e-^ /2 sm9 p e x ^ e ^4>y^ e ~^ y /2 e tk RV 
E y - =E xy (V* /2 cos0 p e z + e" 4 ^ 2 sin^e,) e -<**«/a e -«^/2 e -<fc*w 
E z =^|( e;c +e,)e^. 

In this complicated set of fields, tfi describes the ellipticity of the lasers traveling in 
the e x -e y plane, the major axes of which are tipped by an angle 9 p from vertical. 
When (f> = 7r/2 all four beams are right-hand circular polarized. 5<p x and S(j)y describe 
relative phase differences between the forward- and counter-propagating beams along 
e x and e y respectively; and lastly, 8(j) xy is an overall phase difference between the 
beams traveling along e x and those traveling along e y (a similar phase difference Scj) z 
exists between the e x -e y beams and the e z beam, however, it amounts to simply 
displacing the system along e z ). 

For this set of fields, the scalar light shift [neglecting a u s (4:E 2 y + E 2 ) energy 
offset] of is 

U(v) = U± [cos(2x') + cos(2y')] + U\\ cos a;' cosy', (36) 

where we have introduced the scalar energies U± = 2u s E xy cos (26 p ) and U\\ = 
8u s E xy cos 2 9 P cos(2iy9_), with ip± = (5(f> xy ± <f))/2. The RWA effective magnetic field 
becomes 

fZ = fl± [cos (a/) sin (z' — ip_) + cos (y f ) sin (z' + (f-)] e x 
+ [cos (a;') cos (z' — ip_) + cos (y') cos (z' + <p-)] e y 



On 



sin(x') sin(y') + 5 e z . (37) 



We defined Q± = 2u v E xy E z cos 9 p /H and f2|| = iu v E xy sin(2</? + ) sin 2 9 p /h and 
introduced a dimensionless detuning 6 = 5/KQu. (Here, E z is linearly polarized, 
so it does not have any contributions to 5, as would be the case for a circularly 
polarized beam [48].) In these expressions, we made the simplifying replacements 
x' = IcrX — 5<f> x /2, y' — ftRy — S<j> y /2, and z' = k^z + n/4. These show that the 
phase-differences between beams traveling along e x and e y give rise only to effective 
spatial displacements leaving the topology of the lattice unchanged; in contrast, the 
phase difference between the e x and e y lasers 5<p xy fundamentally change the coupling. 

Somewhat more subtly, transforming the complete Hamiltonian according to the 
unitary rotation U{x)HU^ (x), with U(x) = exp [iF z z' /hj completely eliminates the 
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z-dependence from the Hamiltonian, but as in Sec. 5.1.1 introduces a gauge field k]^F z 
for motion along e z . Under this transformation the effective Zeeman term becomes 

u[x) [f ■ n(z')] u\x) = f • rt{z' = o). 

Therefore the Hamiltonian separates into a sum of independent contributions for 
motion along e z and motion in the e^-e^ plane; without loss of generality, we take 
z' = 7r/4. The expression for fi then reduces to that of Eq. (25) for the physical 
parameters <j) = ir/2 (circularly polarized beams in the e x -e y plane), 
5 = 0. 

With the replacement (3 = f2i|/f2j_, the resulting adiabatic orbital field is 



= 0, and 



P (f*fy~ l~hx9y) sin (-2<p-) 

■ -| 3/2 ' 

/ x 2 + /2 + 2f x f y cos(2^_) + 2 (s + g x9 y 2 



(38) 



This implies that practical implementations of flux lattices require active stabilization 
of the phase between beams traveling along e x and e y , but not the e z beam. For the 
choice = — 7r/4 and 5 = this reduces to Eq. (30). 



Given the dependance of J3(r) on so many parameters, we now consider the first 
order sensitivity to perturbations in 6 = A5 and (f- = — it/ 4 + Aip^; since changes 
in phase sum ip + = 7r/4 + A<^ enter into fiy quadratically, they may be neglected 
at first order. Additionally, the polarization angle dp is generally static in the lab, 
and an imperfect setting can be accounted for by changing the intensity of the Raman 
beams. We find the scalar light shift is unchanged, but the effective coupling becomes 



ft = Q± [cos (x') + Aip_ cos (y')] e x 
+ il± [— A<y9_ cos (x') + cos (?/)] e y 

+ Oil sin(x') sin(y') + AS e z . 



(39) 



Given this, it is surprising but delightful, that we arrive at an orbital field which is un- 
altered at first order. We observe that AS usually results from noise in the magnetic 
field; here this noise must be small compared to the coupling strength Q» , not the 
generally much smaller width of the Bloch bands. 



5.2. Band structure and Chern numbers 

The adiabatic arguments show that flux lattices give rise to large orbital magnetic 
fields with non-zero average. As we learned above, the spatial locations with largest 
magnetic field are also associated with a repulsive adiabatic scalar potential W(r), 
suggesting that without a compensating term from the scalar potential U(r) the 
magnetic field might not be important for atoms in the lowest bands. To address 
this question, we studied the resulting band structure and identified when the Bloch 
bands have non-zero Chern number, in analogy with the band structure of charged 
particles in a magnetic field. We directly compute the band structure from potential 



terms in Eq. (33) combined with the contribution from the usual kinetic energy term, 
and then extract the Chern numbers using the prescription in Ref. [69] . 
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Figure 7. Band structure, (a) Band structure for the lowest four bands showing 
the Chern numbers C n , and the modest energy gap Ai?oi between the ground 
and first excited band. These were computed for Q± = 1.905-El, C|| = 5.1-El, 
U±_ = -1.95E L , Un = 0, 5 = and the phase <f>- = — 7r/4. (b) Ground band 
Chern number as a function of Uj_ and f2ii for the same [7j_ and J7n as above. 
The red contours mark AEqi in this plane, showing the modest parameter region 
where it is non-negligible, and the blue cross locates the maximum gap where the 
band structure of (a) was computed. 



The Hamiltonian described by Eq. ( 33 ) has apparent primitive lattice vectors 
Ui = 2ii/k^e x and U2 = 2ir/kne y (each of these is twice larger than usual for a lattice 
formed by retro-reflected lasers). To compute the band structure in the most simple 



manner we first rotated the coordinate system in Eqs. (36 1 and (37) by 7r/4 in the 
e x -e y plane and defined scaled coordinates x" = y' + x' and y" = y' — x' . In analogy 



with the procedure described in Sec. 5.1.1 we applied a spatially dependent rotation 



U (r) = exp 



i (x" + y") F z 

2h 





ikRyF z 


= exp 


h 



which introduced a gauge term in the kinetic energy. In the example given in 
Sec. |5.1.1| this process completely removed the Hamiltonian's spatial dependence; 
here it does not, but the area of the unit cell is halved [the primitive lattice vectors 
Ui = n/2kn(e x + e y ) and 112 = n/2kn(~e x + e y ) expressed in the initial coordinate 
system are reduced in magnitude by a factor of l/y/2], doubling the area of the 
Brillouin zone [48] . The resulting Hamiltonian has contributions 



2m 



(k x - k L F z /2) + (k y - k L F z /2 



U(r) =U X [cos(x" + y") + cos(a;" - y")] + -j- [cos (x") + cos (y")] 

f • n =fsi z + P+n_ + P_n+. 



(40) 



Where F 4 



F x HF y are 



the usual angular momentum raising and lower operators; 
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the coupling expressed in the helicity basis is quite simple with 

- cos(V ) + cos(y") + 5 

e i( x "+V>-) i e i{y"+V-) i gi(x"+j/"-V3-) , e ~ is P- 

All of these expressions fully respect translational symmetry in the reduced unit cell 
whose reciprocal lattice vectors have magnitude kh — V^kji, and a recoil energy 
El = 2Er. From this, computation of the band structure and it's Chern numbers is 
straightforward. 

Figure [7] depicts the outcome of this computation for an optimally chosen 
parameter set (values given in the caption). While the Chern number of the lowest 
band Co is non-zero over a wide range of parameters, but red contours illustrate 
the most significant limitation for practical implementation of these flux lattices: the 
relatively small gap between the ground and first excited band A£oi. For our optimal 
parameter set, we find a maximal gap of just AE i = 0.107i?L = 0.214ER, far less 
than the U ~ lE^ on-site interaction energy in typical 3D optical lattices (a slight 
improvement is possibly by tuning the quadratic Zeeman term which was absent in 
these computations). This implies that interactions will hybridize several of the lowest 
bands in a way that cannot be described as a perturbation of the lowest band as is 
possible in conventional optical lattices where Ai?oi > lOE^. 

6. Concluding remarks 

We have explored the optical flux lattices produced for ultra-cold atoms in the 
radiation field when both the atom-light coupling and the detuning exhibit an 
oscillatory behavior. We have analyzed not only the magnetic flux but also the 
geometric vector potential generating the flux, as well as the accompanying geometric 
scalar potential. We showed how to deal with the gauge-dependent singularities of 
the AB type appearing in the vector potentials for the optical flux lattices. We have 
present a way to calculate the continuous magnetic flux through the elementary cell via 
the singularities of the vector potential inside the cell. The analysis is illustrated with 
a square optical flux lattice. We have presented a way of creating such a lattice using 
the Raman transitions induced by a set of properly chosen polarization-dependent 
standing waves propagating at a right angle and containing a time-phase difference. 
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